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ABSTRACT 



Context. Outflows from asymptotic giant branch (AGB) and red supergiant (RSG) stars inject dust into the interstellar medium. 
The total rate of dust return provides an important constraint to galactic chemical evolution models. However, this requires detailed 
radiative transfer (RT) modeling of individual stars, which becomes impractical for large data sets. An alternative approach is to select 
the best-fit spectral energy distribution (SED) from a grid of dust shell models, allowing for a faster determination of the luminosities 
and mass-loss rates for entire samples. 

Aims. We have developed the Grid of RSG and AGB Models (GRAMS) to measure the mass-loss return from evolved stars. The 
models span the range of stellar, dust shell and grain properties relevant to evolved stars. The GRAMS model database will be 
made available to the scientific community. In this paper we present the carbon-rich AGB model grid and compare our results with 
photometry and spectra of Large Magellanic Cloud (LMC) carbon stars from the SAGE {Surveying the Agents of Galaxy Evolution) 
and SAGE-Spec programs. 

Methods. We generate models for spherically symmetric dust shells using the 2Dust code, with hydrostatic models for the central 
stars. The model photospheres have effective temperatures between 2600 and 4000 K and luminosities from ~2000 L Q to ~40000 L Q . 
Assuming a constant expansion velocity, we explore five values of the inner radius R iu of the dust shell (1.5, 3, 4.5, 7 and 12 /? star ). We 
fix the outer radius at 1000 R[ n . Based on the results from our previous study, we use amorphous carbon dust mixed with 10% silicon 
carbide by mass. The grain size distribution follows a power-law and an exponential falloff at large sizes. The models span twenty-six 
values of 1 1.3 jim optical depth, ranging from 0.001 to 4. For each model, 2Dust calculates the output SED from 0.2 to 200 fim. 
Results. Over 12000 models have dust temperatures below 1800 K. For these, we derive synthetic photometry in optical, near- 
infrared and mid-infrared filters for comparison with available data. We find good agreement with magnitudes and colors observed 
for LMC carbon-rich and extreme AGB star candidates from the SAGE survey, as well as spectroscopically confirmed carbon stars 
from the SAGE-Spec study. Our models reproduce the IRAC colors of most of the extreme AGB star candidates, consistent with the 
expectation that a majority of these enshrouded stars have carbon-rich dust. Finally, we fit the SEDs of some well-studied carbon stars 
and compare the resulting luminosities and mass-loss rates with those from previous studies. 

Key words. Stars: AGB and post-AGB, Radiative transfer, Stars: carbon, Stars: mass-loss, (Stars:) circumstellar matter, (Galaxies:) 
Magellanic Clouds 



1. Introduction 

During the final stages of their evolution, low- and intermediate- 
mass stars (0.8 to 8 M ) ascend the asymptotic giant branch 
(AGB). AGB stars are prominent members of stellar popula- 
tions of ages ~0.2-2 Gyr. Characterized by mass-loss rates of 
up to 10~ 4 M yr -1 , stars in the AGB phase inject a signifi- 
cant fraction of their mass into the interstellar medium (ISM). 
This mass loss is thought to occur in two steps: stellar pul- 
sations first levitate material to the cool, outer layers where 
dust grains form. Interaction with stellar photons then accel- 
erates the dust grains, which in turn d r ag the gas along with 
them (e.g.. iGoldreich & Scovilld 119761: iHofner & Dorfilll997t 
IWachter et al.ll2002t lHofnenl2009l) . AGB stars of masses ~1- 
4 M undergo the third dredge-up process (e.g.. Ilbenl 1 1983: 
iKarakas et alJ l2002h which transports the products of nuclear 



* The model grid is available at the CDS via anony- 
mous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http : //cdsweb . u-strasbg . fr/cgi-bin/qcat? J/A+A/ 



reactions, including carbon, into the outer layers. Each dredge- 
up event increases the carbon abundance relative to oxygen un- 
til eventually the C/O ratio exceeds unity and carbon stars are 
born. Car bon dust is more efficient a t absorbing optical pho- 
tons (e.g.. IWallerstein & Knappl ll998) and it has higher emis- 
sivity at infrared (IR) wavelengths. High rates of mass loss from 
C-rich AGB stars make them majo r contributors of atomic car- 
bon and carbonaceous dust g rains dDwekl [T998: Matsuura et al.l 
2009; ISrinivasan et alj |2009) to the ISM and may eventually be 
assimilated into star-forming regions. The thousands of carbon 
stars present in galaxies with intermediate-age stellar popula- 
tions contribute substantially to their integrated bolometric an d 
near-IR (NIR) luminosities dFrogel et al.ll 1990b iMarastonfl 998 ), 
and this contributio n is somewh at higher at lower metallicity 
(e.g., Fig. 13 in iMarastonl |20 05 ) . Therefore, we must quantify 
the carbon-star dust output in order to study the dust cycle in 
galaxies as well as constrain stellar population synthesis mod- 
els. 
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The study of Milky Way AGB stars is inhibited by the 
presence of substantial interstellar extinction, and large un- 
certainties in distance determinations. The Large Magellanic 
Cloud (LMC) offers the combination of proximity (~50 kpc; 
Ivan Leeuwen et all 2007b. low line-of- sight extinction (E(B - 
V) ~ 0.07 5 mag; Schlegel et al. 1998) and favorable orienta- 
tion (~24°; Weinberg & Nikolaev 2001). These properties allow 
in-depth studies of the entire LMC AGB population. O ne such 
study, the Spitzer Space Telescope (I Werner et al.ll2004l) Legacy 
program SAGE (Sur veying the Agents of a Galaxy's Evolution; 
iMeixner et al. 2006), imaged a 7x7° area centered on the LMC 
and found over 6 million p oint sources, inc l uding thousands 
of ca rbon-star candidates (Blu m et all 120061: [Srinivasan et alJ 
2009). Follow-up spectroscopy from the Infrared Spectrometer 
(IRS. lHouck et al.ll2.QQ4j) on Spitzer was o btained as part of the 
SAGE-Spec program dKemper et al.ll2010l) . 

SAGE provides an ideal dataset for AGB studies. 
ISrinivasan et al.l (12009) (hereafter, Paper I) calculated mid-IR 
(MIR) excess fluxes for AGB candidates identified from SAGE 
photometry and used these to estimate the total dust injection 
rate into the LMC. A more precise estimate for the injection rate 
requires detailed radiative transfer (RT) modeling of the spec- 
tral energy distribution (SED) of each star in the candidate list. 
Many authors have c omputed such detailed models for LMC 
carbon stars (see, e.g., van Loon et alJ[T999l iGroenewegen et ail 
2009; ISrinivasan et all 120101) . The computation of individual 
models becomes time-consuming for large samples such as the 
SAGE dataset. For this purpose, an alternative would be to 
compare the observed SEDs to those of pre-computed models. 
Ideally, such a grid of models should account for the photo- 
spheric absorption due to atomic and molecul ar species (e.g., 
Gautschv-Loidl et al 1 12004b lAringer et al.ll2009h and dynamical 
effects such as pulsation-dr iven shocks (e.g., iHofneret al.ll200H 
Matt sson et al.ll2007l l2010h. as well as the repr ocessing of stel- 
lar radiation by dust (e. g.. lGroenewegenll20 06). The effects of 
stellar evolution can be folded into the grid either by performing 
the above calcul ations on stars sampl ed along AGB evolutionary 
tracks (see, e.g jMattsson et alj2007h or by folding the effects of 
dust into stellar population mo deling (e.g.. iMarigo et al. 2008; 
iGonzalez-Lopezlira et alJl20loh . 

The method of comparison with model grids is suitable for 
quickly constraining the general properties of large photometric 
datasets such as the LMC SAGE set. With the intention of quick 
SED fitting in mind, we have generated the Grid of RSG and 
AGB ModelS (GRAMS). The GRAMS grid consists of radia- 
tive transfer models computed using model stellar photosphere 
fluxes and assuming various values of stellar and dust shell pa- 
rame ters. While some m odel grids are available in the literature 
(e.g.. lGroenewegenll2006l) . they typically use only a few tem- 
plate solar-metallicity photospheric spectra as the basis for the 
grid. Our models cover a large region of the parameter space 
spanned by AGB stars. One of the motivations for our study 
is the recent availability of large sets of photosphere models 
with an improved tr eatment of molecular spectral features of 
low- metallicity stars dGautschv-Loidl etal ] l2004t |Aringer et al. 
I2009h . This enables us to probe the dependence of mass loss over 
a large range of stellar parameters. Other studies typically equate 
the dust temperature at the inner radius to the dust condensation 
temperatu re, which is usually fixed at 1000 K for ca rbon dust 
(see, e.g., Ivan Loon et al.l [19991; IGroenewegen! |200oT) . We cir- 
cumvent this assumption by adopting a more general treatment 
- we specify the inner radius as input to the modeling and this 
automatically determines the temperature of the dust in the shell 
when the model is computed. The output SEDs are very sensitive 



to the inner radius R m (see, e.g.. ISrinivasan et aT1l2010l) and we 
incorporate this dependence in our grid by computing models for 
different Ri n values. We thus provide a large grid of models that 
is complementary to the currently available grids constructed by 
other authors. In previous pape rs in this series, we d etermined 
the properties of oxygen-rich (Sarge nt et al. 1 1201 Oal her eafter 
Paper II) as well as carbonaceous (ISrinivasan et alJl2010l here- 
after P aper III) dust grains for use with the grid. IGroe newegen 
(2006) incorporate different dust species in their grid. While the 
current incarnation of our grid is computed for a fixed set of dust 
properties, we will investigate other types of dust in future ver- 
sions. 

Fitting the SEDs of our entire dataset will allow us to inves- 
tigate the LMC AGB mass-loss return. In this paper, we provide 
the details of the GRAMS carbon -star grid. Our O-r ich grid is 
described in a companion paper (Sarg ent et alJl20lH hereafter 
Paper IV). The remainder of the paper is organized as follows. 
We detail the parameter selection for the carbon-star grid in 
Sect. |2] and describe the computational procedure and the cal- 
culation of synthetic photometry for the resulting grid in Sect. [3] 
In Sect. [4] we compare the GRAMS synthetic photometry with 
SAGE observations on color-magnitude as well as color-color 
diagrams. As a further validation of the grid, we perform sim- 
ple chi-squared fits to the SEDs of some sources that have been 
previously studied in detail. We present our conclusions in Sect, 
i 



2. The Model Grid 



We use the 2Dust code dUeta & M eixner 2003) to populate the 
grid. 2Dust is a radiative transfer code for axisymmetric sys- 
tems. While the mass loss is spherically symmetric for most of 
the AGB phase, the highly evolved AGB stars as w ell as their 
post- AGB successors exh ibit bipolar geometries (e.g JUeta et al.l 
2000; IMeixner et alj|2002h . The dust around post- AGB stars is 
also similar to AGB star circumstellar dust. We would like to 
produce models for these objects as well, but in this paper we 
assume spherically symmetric shells for simplicity. We will con- 
sider the effect of non-spherical geometries in future versions of 
the grid. In order to solve the radiative transfer equation, 2Dust 
discretizes the dust shell into a 2-dimensional polar grid. The 
code then computes the radiation field at each grid point, dis- 
cretizing the field into a set of incoming rays ("characteristics") 
that converge on the grid point from all directions. During each 
iteration, the code ensures self-consistency by requiring a global 
luminosity constancy throughout the dust shell. For a given stel- 
lar spectrum, dust shell geometry (inner and outer radius, den- 
sity variation) and a set of dust grain properties (species, optical 
depth at a reference wavelength, grain size distribution), 2Dust 
calculates the luminosity^] and mass-loss rate for the system. The 
code then solves the radiative transfer for the output star+dust 
spectrum. In this section, we discuss our parameter selection 
and provide details of our computational procedure. Table[T]lists 
some of the input/output parameters as well as the range in each 
parameter covered by the GRAMS carbon-star grid. 



1 Calculation of the luminosity requires the distance to the star, which 
we set at 50 kpc. 
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Table 1. GRAMS carbon-star model grid parameter coverage 



Photosphere model" 

L«(L ) -1100 to -26000 

T eff (K) 260010 4000(100)* 

loggtcms" 1 ] -1.0 to 0.0 (0.1 ) e 
M(M ) 1,2, 3 and 5 

C/O 1.4, 2.0 and 5.0 

Dust shell properties 

Ria(R,) 1.5,3,4.5,7,12 
Ret (Rin) 1000 
density profile p(r) oc r~ 2 
v exp (kms _1 ) 10 

Dust grain properties 

Species AmC d +SiC e 

SiC fraction 10% 

t(1 1 fim) 10~ 3 to 10" 1 (5 per dex), 

0.1 to 1 (0.1) and 

1.5 to 4 (0.5) 

Size distribution KMH / 

flminOum) = 0.01 
fl (/"m) = 1 
r = 3.5 

Mass-loss rate and dust temperature 

MdusttMoyr- 1 ) 1.5 x 10" 12 to 2.1 x 10~ 7 
M gas (M yr- 1 )* 3.0 x lO" 10 to 4.3 x 10~ 5 



'gas 

T in (K) 710 to 1800'' 

Notes. <D) Photosphere model from Arin ger et alj ([2009) (A09). <6> Where applicable, parameter increments are supplied i n parentheses . M N o 
A09 model photospheres were available for logg = -0.1 (rf) Am orphou s carbon grains, p = 1.8 g cirT 3 , optic al constants from lZubko et all fl 996). 
w SiC grains, p = 3.22 g cm 4 , optical constants from Pegourie ( 1988). Size distribution from Kim et al. ( 1994): 

n(a) ~ aT y exp (-a/ao) with a > a m j„. te) Assuming a gas:dust ratio of 200. (/,) As explained in the text, the grid consists of models with temperatures 
cooler than 1800 K. 



2. 1 . Input for the grid 
2.1.1. Stellar spectrum 

We represent the central stars using the lAringer et alJ ([2009) 
COMARCS hydrostatic models of AGB star photospheres. 
Aringer et al. (hereafter, A09) take into account the contribution 
to the emergent spectrum from atomic and molecular absorption 
lines, which can cause the optical and near-infrared spectrum 
to deviate significantly from that expected from a blackbody at 
the same effective temperature. The A09 set consists of spheri- 
cally symmetric COMARCS hydrostatic models calculated with 
opacities generated for CO, C2, CN, C2H2, HCN, C3 and a num- 
ber of other molecules using the COMA (Copenhagen Opacities 
for Model Atmospheres) code. The parameter coverage of the 
models reflects the range of effective temperatures, surface grav- 
ities and C/O ratios predicted by synthetic evolution models. 
These are the most accurate models available at present. A09 
find good agreement between their models and data for warmer 
AGB stars. They do not account for dynamical effects and the 
processing of starlight by dust. In addition to an entire grid calcu- 
lated for solar metallicity, there are subsets of models for metal- 
licities of Z=0.33 and 0.1 Z , corresponding to th e Magellanic 
Clouds (Zt,mc/Zp)~0.3-0.5 and Zsmc /Zn~0. 1-0.2: iDufour et alJ 
[1983 [Bernard etalJl2008h . As we are interested in first repro- 
ducing the observed range of colors for the LMC, we selected 



only their Z = 0.33 Z subset of 131 models of masses < 5 M . 
We briefly outline the range of values covered by this subset. For 
details on the entire grid, we refer the reader to Sec 2.1 of A09. 

The effective temperatures of the Z=0.33 models (hereafter, 
the LMC set) range from 2600 K to 4000 K in increments 
of 100 K. While carbon stars may have lower effective tem- 
peratures than this range covers dMattsson et alJ 120081) . the 
A09 mod els provide adequate co verage for the LMC C-rich 
sample (Groenewegen et al. 2009]). As noted in A09, most 
models with temperatures hotter than 3500 K will not be on 
the AGB, but they are included in order to model carbon-rich 
objects in the post-AGB phase. The surface gravities range from 
log (g[cm s -1 ]) = to -1 in steps of 0.1, but not all logg values 
are available for a given temperature. In general, the cooler 
the temperature, the lower the minimum logg available. The 
LMC set was calculated completely for a stellar mass of 2 M , 
resulting in 126 models. There are also three 1 M models with 
r eff =2600 K and 3000 K, as well as one r eff =3000 K model 
each for mass of 3 M and 5 M . A09 find that a change in 
stellar mass has only a minor effect on the near-IR color and 
bolometric magnitude (see their Fig. 8). These models have 
C/O ratios of 1.4, 2 and 5. Higher C/O ratios are expected 
at lower Z due to the under-abundance of oxygen. Due to 
this fact, the C/O ratio can become significantly higher than 
unity within the first few dredge-up events. For the range of 
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model parameters corresponding to each temperature range, 
see Table 1 in A09. The luminosities of the selected models 
are roughly in the range 1 100 - 26 000 L (Af bo i = -2.9 to -6.3). 

In Fig. [2] we compare this range to the luminosity distribu- 
tion calculated for S AGE C-rich and extreme AGB candidates. 
Mari go et al.l Jl999) estimate that the observed LMC carbon-star 
luminosity function extends between M\, \ = -3 and -6.5; these 
numbers are in good agreement with ourfl IRAS 04496-6958 
is the most luminous spectroscopi cally confirmed ca rbon star in 
the LMC (> 30000 L ; see, e.g. ISpeck etal]|2006h . Based on 
SAGE photometry alone, we estimate its luminosity to be about 
35 000 L e (Mboi ~ -6.7) which is still below the classical lumi- 
nosity limit. Fig.|2]also shows the GRAMS luminosity coverage. 
To fit the SEDs of the most luminous carbon stars, we may have 
to scale up the luminosities of our models in a procedure similar 
to that adopted in Sargent et al. (Paper IV) for RSGs and O-rich 
AGB stars. Our aim in constructing the GRAMS model grid is 
to span the range of parameters observed for LMC carbon stars; 
it is possible that in this attempt we may generate some mod- 
els with extreme luminosities that are not representative of real 
carbon stars and may never be matched to data. We do not use 
luminosity scaling in the current version of the carbon star grid; 
however the scaling is easy to apply once the model grid has 
been generated, and we will consider such a scheme when we 
perform SED fitting of the SAGE data in upcoming papers. 

2.1.2. Shell geometry 

We assume a constant mass-loss rate at a constant outflow ve- 
locity u exp , leading to an inverse-square de nsity distribut ion in 
the shell. Models of stationary winds (e. g. . IWoit ke 2006) show 
that the wind is accelerated from rest near the stellar surface to 
an almost constant speed within a few stellar radii. The accel- 
eration of the wind therefore affects the density profile in the 
regions of the dust shell that contribute significantly to the opti- 
cal extinction and mid-IR emission. In the current paper, we do 
not take this radial dependence of the outflow velocity into ac- 
count; ignoring such complications enables us to generate a full 
grid spanning the observed range for a small set of parameters 
while allowing us to compensate for the simplifications at a later 
time. We also ignore any dependence of the expansion velocity 
on the metallicity or gas:dust ratio. As the dust mass-loss rate 
calculated by the code is directly proportionate to the value of 
u exp chosen (see, e.g., Eq. 2 in Paper III), it is straightforward 
to incorporate the effect of changing u exp once the grid has been 
populated. Following our discussion in Paper III, we use u exp = 
10 km s _! for all our models. 2Dust accepts a user-defined den- 
sity function, which can be used to study time-dependent mass 
loss. We will also consider this possibility in future versions of 
the grid. 

As discussed in Paper III, the mass-loss rate and output 
spectrum are very sensitive to the value of the inner radius 
R m . The inner radius determines the hottest temperature of the 
dust, which is one of the 2Dust output parameters. We calcu- 
late models for Ri n = 1.5, 3, 4.5, 7 and 12 stellar radii. Simple 
energy-balance estimates suggest that amorphous ca rbon dust 
should form within a few stellar radii (iHofnerl 12007). A lower 
limit of about 1.3 R in is suggested for SiC formation in ex- 



2 We calcu late luminosities higher than the classical AGB limit 
(M bo i = -7.1, Paczvnski 1970) for a handful of our extreme AGB can- 
didates; while a majority of these are well-studied O-rich AGB stars, 
there are no carbon stars in this subsample. 



treme carbon stars bv lSpeck et al.l (120091) . The range considered 
here also agre es with observation s of Galactic carbon stars (e.g., 
IRC+10216; iDanchiet all 119951) and result s from RT model- 
ing of LMC stars ( e.g.. Ivan Loon et al1ll999l) . We note here that 
2Dust uses the inner radius as input to determine the dust tem- 
perature (Ti n ) as a function of radius in the shell. This means 
that we cannot directly restrict the range of T[ n before the grid is 
computed. However, we filter out models with high and/or un- 
physical dust temperatures from the grid once it is generated (see 
Sect. 1531 for details). 

The outer radius determines the total amount of mass in the 
shell and the mass-loss timescale. While these are important 
quantities, we focus on obtaining mass-loss rates, which are only 
weakly sensitive to changes in the value of the outer radius. We 
ensure that the outer radius is large enough so that we do not 
miss any contributions to the flux from the outermost regions. 
While modeling the shell around OGLE LMC LPV 28579 in 
Paper III we found that a value of R out = 1000 Ri n was suffi- 
ciently large to satisfy this condition. We use the same value for 
all the models in our grid. For u exp =10 km s and R out = 1000 
R m , the mass-loss timescales corresponding to the smallest and 
largest circumstellar envelopes in our grid are 548 yr and 19 620 
yr respectively. 

2.1.3. Dust grain properties 

We modeled the SAGE photometry and SAGE-Spec spectrum 
of OGLE LMC LPV 28579 in Paper III for the purpose of se- 
lecting a set of carbonaceous dust grain properties for GRAMS. 
As a result, we choos e a mixture of am orphous carbon (AmC; 
optical constants from Zubk oet al Jl 1 996|) with 1 0% silicon car- 
bide by mass (optical constants from Pegourie 1988). As input, 
2Dust also requires the optical depth specified at a reference 
wavelength. In Paper EH, we used the SiC feature observed in 
the SAGE-Spec spectrum of LPV 28579 to constrain the opti- 
cal depth as well as the SiC content of our dust model, so it 
was convenient to specify the optical depth at 11.3 /mi. We fol- 
low the same practice in this paper, because the dust composi- 
tion used in the grid is identical to that of Paper IIlQ We will 
revise the current convention when we incorporate more dust 
types into our grid. The values of rn.3 in our grid range from 
10~ 3 to 4. We consider five optical depths per decade between 
10~ 3 and 0.1 (1, 2, 4, 6 and 8 xl0~ 3 and so on), and ten values 
between 0. 1 and 1 . Additionally, we calculate models for Tn.3 = 
1.5, 2, 2.5, 3, 3.5 and 4, bringing the total to twenty-six unique 
optical depth values. The grain sizes ar e distributed acco rding 
to n(a) oc cT y e~ alaa (KMH distribution. iKim et all 1 1994 with 
power-law index y = 3.5, minimum grain size a m i a = 0.01 /im 
and exponential scale height oq = 1 /an. For the AmC: SiC mix- 
ture considered here, the average grain siz^|is ab out 0. 1 /an, the 
value typic ally used in s i ngle-s ize models (e.g., iGroenewegenl 
120061) . The iMafhis et al.1 (1 19771) (MRN) distribution, given by 
n(a) oc aT y for a m j n < a < a max places a strict limit on the 
maximum grain size. In this sense, the KMH model represents 
a more realistic grain size prescription. In Paper III, we showed 
that the output spectrum showed only a weak dependence on the 
exponential factor «o of the KMH model beyond ao ~ 0.1/an. 
Our choice of the KMH model thus means there is one less 



3 For reference, the optical depth at 1 fim is about 10.64 times its 
value at 1 1.3 //m for the chosen set of dust properties. 

4 For the GRAMS grid, averages over grain size space are computed 
by wei ghting according to gra in surface area (Harrington averaging 
scheme, Harring ton et al.lll988h 
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parameter to constrain. Owing to its maximum grain size re- 
quirement, the MRN distribution has a smaller average grain 
size. This results in higher absorption and therefore lower opti- 
cal flux and more mid-I R emission. This effect was discussed in 
lUeta&Meixnerl (12003). In our study, we computed 2Dust mod- 
els for a few GRAMS models using the MRN prescription and 
found that the maximum difference in fluxes is less than 5%. 



3. Generating The Grid 

3.1. Preparation 

We use 131 A09 photospheres and we explore five values for the 
inner radius as well as twenty-six optical depths, making a total 
of 17 030 output models. Generating such a huge grid of models 
requires that 2Dust be run in the non-interactive mode, with all 
the required information provided in the form of various input 
files (For details of the input format, we refer the user to Sect. 
3.1 of the 2Dust manual Q) 

The A09 photospheric spectra comprise of fluxes at about 
8 000 wavelengths ranging from ~0. 44-25 fim. As we are in- 
terested in providing synthetic photometry over all the broad- 
band filters available to us (optical U through MIPS24), we ex- 
trapolate the A09 fluxes onto a larger wavelength grid. On the 
long-wavelength side, we extend the grid beyond the MIPS24 
band, to allow for future comparisons with MIPS 70 and 160 fim 
fluxes or other long-wavelength data, if available. However, this 
raw extrapolated spectrum cannot directly be fed into 2Dust; the 
computational speed of each run is approximately linearly pro- 
portional to the number of points in the wavelength grid (Sect. 
3.1.5 of the 2Dust manual). We therefore sample the extrapo- 
lated spectrum at about 130 wavelengths ranging from 0.2 fim 
to about 200 fim. We sample more near- and mid-infrared wave- 
length points as we desire to fit the mid-IR photometry and spec- 
tra available as part of the SAGE and SAGE-Spec programs; the 
constraint on the total number of wavelength points then reduces 
the number of samples for A < 1 fim. As a result of this under- 
sampling, we expect that we may not reproduce many of the 
sharp atomic/molecular features observed in the optical spectra 
of AGB photospheres, which would affect the heating of dust 
grains as well as the broadband optical colors of the models 
with optically thin dust shells. This latter issue is of less con- 
cern since we are interested in reproducing the near- and mid-IR 
photometry. However, in order to have a more accurate treat- 
ment of dust heating, we will incorporate a better sampling in 
the optical in future versions of the grid. Fig.[T]demonstrates the 
effect of undersampling on the output model for a r e jf = 3300 
K, log(g[cm s -1 ]) = -0.2 photosphere with an optically thin 
(th,3 = 2 x 10 3 ) dust shell. While the resulting GRAMS model 
does not reproduce all the narrow optical features, it shows some 
of the large-scale variation in the spectrum. The sampling de- 
scribed above also changes the integrated flux calculated for 
each model by at most a few percent - we found that the lumi- 
nosities calculated from the raw and sampled spectra agreed to 
within 5% for 115 of the 131 A09 photospheres, with the max- 
imum discrepancy being around 8%. For consistency, we only 
use the post-sampled luminosities to tag our resulting models. 

The number of discrete radial zones in the shell (NRAD) di- 
rectly determines the total memory allocated to the code. Based 
on our tests, we chose to have 64 radial zones for all of our mod- 
els. The code converges to a solution when, at each grid point, 



5 The 2Dust manual is available at 
http://www. stsci . edu/science/2dust/2dust_manual .pdf . gz 



the fractional change in the integrated flux between consecutive 
iterations is less than a user-defined tolerance level. This level is 
specified in the code in the form of the CONDITION parame- 
ter, which we fixed at 5xl0 -4 for all the models. The execution 
time is then most sensitive to the optical depth. There is a trade- 
off between output precision and execution efficiency, which are 
regulated internally by two parameters: one (VSPACE) sets the 
smallest line integration step size in terms of the local mean free 
path length, while the other (MXSTEP) sets the maximum num- 
ber of allowed steps for line integration along each characteris- 
tic. The smallest step size is inversely proportionate to VSPACE 
value chosen, so large VSPACE values are required for higher 
precision in the results. However, VSPACE determines the actual 
number of steps along the characteristic, which cannot exceed 
MXSTEP. This latter value controls the array length for long 
characteristics and therefore the memory access time upon exe- 
cution of the code. Therefore, we need to optimize the choices 
for MXSTEP and VSPACE as a function of optical depth in or- 
der for the code to converge in a reasonable amount of time. For 
this purpose, we ran a number of test models on a MacBook Pro 
laptop with a 3.06 GHz Intel Core 2 Duo processor and 8 GB 
RAM. The test models sampled optical depths over the entire 
range of values considered for the grid. These models explored 
various (MXSTEP, VSPACE) pairs, enabling us to select ideal 
values over different optical depth ranges. This information is 
incorporated into the input files for each model in the grid, and 
this input enables us to run 2Dust in the non-interactive mode in 
an automated fashion. 

2Dust tracks the temperature of the dust grains during each 
iteration. Assuming radiative equilibrium, the code uses these 
temperatures to obtain the corresponding intensities k v B v from a 
look-up table, which is computed on execution. This table con- 
sists of a temperature grid whose lower and upper limits are 
TBOT = 2.7 K and TTOP = 1000 K (the typical value chosen 
for the condensation temperature of carbon dust) by default. It is 
not easy, however, to constrain r con d as it depends on many fac- 
tors such as the C /O ratio and the gas pressure {e.g., Sect. 3.4 of 
Speck et al. 2009, and references therein). A detailed treatmen t 
of r con d is beyond the scope of our work. [Speck et al. (2009) 
suggest that carbon dust can form at temperatures above 1400 K 
even at low mass-loss rates. They also find that for high C/O ra- 
tios, graphite grains can form at temperatures of ~1800 K. While 
our current models only consider amorphous carbon dust grains, 
we are interested in as few constraints on the output parameters 
as possible so as to be able to extend our treatment by including 
more grain types in the future if required. Therefore, we assume 
an upper limit of 1800 K for r con( j. We allow for slightly higher 
dust temperatures in the look-up table by setting TTOP to 2000 
K to avoid convergence issues (see Sect. 13 . 3b , and we enforce 
the 1800 K constraint on r conc j once the grid is generated. 

3.2. Batch job submission 

We automated the 2Dust execution using the Magique su- 
percomputing cluster at the Institut d'Astrophysique de Paris. 
Magique is a cluster of 96 AMD operton nodes interconnected 
with a Gigabit ethernet network, 84 of which are reserved for 
scientific computation. Of these, we used the bisockets single 
and double core nodes (18 and 64 nodes respectively), which al- 
lowed 8 GB RAM per job. 2Dust was compiled using the Intel 
FORTRAN compiler ifort available on Magique. The batch 
queue system on the cluster is managed by PBSpro v8.0, which 
automatically allocates a position in the queue for submitted jobs 
and continuously monitors their status. Our jobs were grouped 
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Fig. 1. Aringe r~et al.1 J20 09) photosphere spectrum (gray) with T e ff = 3300 K and log (g[cm s -1 ]) = -0.2, showing many strong 
features at shorter wavelengths. The vertical dashes show the wavelengths at which this spectrum was sampled for input to 2Dust. 
The resulting GRAMS output model with optical depth T\\ .3 = 2 x 10~ 3 (black) incorporates some of the variation in the optical. 



according to the A09 photosphere used and they were submit- 
ted for execution in batches consisting of up to five sets of pho- 
tospheres (i.e., upto 650 models), depending on the queue sta- 
tus and the load on the cluster from other users. Depending on 
the optical depth, the time taken for a single model to converge 
ranged from under a minute to about two hours. We were able to 
generate the entire grid in under two weeks. 



3.3. Filtering out models based on dust temperature 



Not all combinations of input parameters detailed in Tab. Q~]will 
result in realistic or even physical representations of the dust 
around carbon stars. For instance, we expect that the output 
dust temperatures for high T^, low and high T11.3 may be 
higher than the dust condensation temperature r con d. If the lo- 
cal radiation intensity is high enough to warm the dust grains to 
temperatures beyond TTOP (the maximum value in the lookup 
table, which was set to 2000 K as described in Sect. 13. U . the 
code is unable to track these temperatures and intensities, re- 
sulting in unphysical (extremely large or negative) dust temper- 
atures in the inner regions of the shell. We found that 1860 mod- 
els did not converge because of this issue. As expected, these 
models consisted mainly of shells that were too close to the star 
(Ri n = 1 .5/? st ar), but also included some models with larger R m 
and T e ff values warmer than 3500 K. The number of converged 
models in our grid is therefore 15 170. We then impose the con- 
straint that the dust temperatures be cooler than the chosen con- 
densation temperature (r con d = 1800 K, see Sect. 13. Il l, which 
eliminates 2915 more models. Our final grid thus consists of 
12 243 models with dust temperatures cooler than 1800 K. The 
number will reduce further if we adopt a lower r con d value; for 
instance, there are only ~9000 models with T m < 1400 K. 



3.4. Publicly available synthetic photometry 

For comparison with photometric data, we convolved the SED 
output from each 2Dust model with the relative spectral response 
(RSR) curves for the various broadband filters. In particular, we 
produced synthetic photometry in this fashion for the MCPS 
UBVI0, 2MASS JHKJ3 and Spitzer IRAC0 and MIPS2/Qbands. 
We compare the synthetic photometry against SAGE data in the 
next section. 

Synthetic photome try was als o derived for the AKARI 
dMurakamiet al.1 120071) and WISE dWrightll2009l) passbandfQ 



6 The MCPS magnitudes were placed on the Johnson-Kron- 
Cousins UBVI system. The detector quantum efficiency (QE) 
curve available on the Las Campanas Observatory website at 



http://www. lco. cl/lco/telescopes-information/irenee-du-pont/i 
was extrapolated to 0.3 pm. The QE was also assumed to drop linearly 
to zero from its value at 0.84 pm, the last wavelength provided on 
the graph, to 1.13 pm, the wavelength corresponding to a photon 
energy of 1.1 eV, which is the energy of the band gap of sili- 
con. The B- and V-band transmission profiles were obtained from 
http://www. lco. cl/lco/telescopes-information/irenee-du 
-pont/instruments/website/direct-ccd-manuals/direct-ccd 
-manuals/3x3-filters-for-ccd- imaging (Harris B profile LC- 
3013 and Harris V profile LC-3009 respectively). The U- and 
I-band filter profiles were obtained from I. Thompson (2009, private 
communication). 

7 The 2MASS filte r relat ive spectral responses (RSRs) de- 
rived by ICohen et alJ d2003h were obtained from the 2MASS 
All-Sky Data Release Explanatory Supplement, available at 
http : //www. ipac . caltech . edu/2mass/releases/allsky/doc/ 
sec6_4a.html . 

8 RSRs available at http://ssc . spitzer.caltech.edu/irac/ 
calibrationfiles/spectralresponse 

9 RSR available at http : //ssc . spitzer . caltech. edu/mips/ 
calibrationfiles/spectralresponse/ 

10 Filter response curves available at 



http : //www . ir . isas . jaxa . jp/ASTRO-F/Observation/RSRF/IRC_FAD/i: 
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The AKARI Sll band (centered around 11 /urn) as well as 
the WISE W3 (centered around 12 pun) band directly sample 
the 11.3 fim SiC feature in carbon stars and the ~ 10 yum sil- 
icate feature in O-rich AGB stars. Ilta et alJ (|2008) presented 
the AKARI survey of the LMC and discussed color-color and 
color-magnitude diagrams. The increase in SiC feature strength, 
therefore, is reflected in the [SI 1]-[L15] color. 

Spectra and synthetic photometry for the GRAMS O-rich 
and C-rich models will soon be available on the 2Dust website 
at the Space Telescope Science InstitutJH. The format of the 
publicly-available data is still undergoing testing and will be re- 
fined through input from the end-users. In addition to the filters 
mentioned above, we can make synthetic photometry available 
for other filters on request. 

4. Results 

In this section, we demonstrate the coverage of various parame- 
ters by the GRAMS grid. We discuss the luminosities and dust 
mass-loss rates. We also compare the GRAMS synthetic pho- 
tometry with data for the SAGE C-rich and extreme AGB can- 
didates as well as SAGE-Spec carbon stars on color-magnitude 
diagrams (CMDs) and color-color diagrams (CCDs). 

4.1. Luminosities and dust mass-loss rates 

The large number of A09 models provide us with a relatively 
dense coverage of the luminosities expected for LMC AGB stars. 
This fact is evident from Fig. [2] which compares the GRAMS 
model luminosities with the distribution of luminosities calcu- 
lated in Paper I from SAGE photometry for C-rich and extreme 
AGB candidate^- The GRAMS models provide excellent cov- 
erage over the entire range of luminosities, except at the faint 
and bright ends. There is a gap visible between the two lowest 
available luminosities of 1100 L and 2200 L . There are 50 
candidates (8 C-rich, 42 extreme) brighter than the brightest C- 
rich GRAMS model available, and 60 candidates (45 C-rich, 15 
extreme) fainter than ~ 2200 L . 

The dust mass-loss rates in the GRAMS grid range from 
1.5x 10~ 12 M yr~' to 2.1 x 10~ 7 M yr _1 . As mentioned in Sect. 
12.1.31 2Dust can only constrain the ratio Mdust/^exp, so the range 
of mass-loss rates produced by the grid depends on the value of 
u exp chosen. The lowest values of mass-loss rates calculated from 
the IR excesses for SAGE carbon star candidates in Paper I were 
around 10~ n M yr _1 . As we are interested in covering as much 
of the parameter space as possible, our grid contains mass-loss 
rates well below this limit. Carbon stars can have dust mass-loss 
rates up to a few times 10 7 M yr~' (e.g.. iGroenewegen et al.l 
2009). Our grid does not currently have adequate coverage of the 
highest luminosities expected for carbon stars (see Sect. I2.1.U . 
but the models can be scaled to these higher luminosities as in 
Paper IV. This would also result in higher mass-loss rates. As 
an example, if we scale the brightest model (L = 2.6 x 10 4 
L ) to the classical AGB limit (L - 5.3 x 10 4 L ) while keep- 
ing all the other parameters constant, the corresponding highest 



and 

http : / /www . astro . ucla . edu/ ~wright/WISE/passbands . html 
respectively 

11 http : //www . stsci . edu/ science/2dust/grams_models . cgi 

12 The luminosities in Paper I (Fig. 4), obtained using the trapezoidal 
rule, were overestimated for the brightest stars because of erroneous 
coefficients in the sum. In the present paper we have used the correct 
coefficients, leading to lower luminosities as shown in Fig. [2] 



mass-loss rate is 3 x 10~ 7 M yr . The total (gas+dust) rate 
depends on the gas:dust ratio, x ¥. As AGB stars produce their 
own carbon, it is easier at lower metallicities to produce an ex- 
cess abundance of carbon relative to oxygen; however, whether 
the gas:dust ratio shows a dependence on meta llicity remains 
an open question (see, e.g.. Ivan Loon et al] |2008, and references 
therein). Using *F = 200, the value determined for Galactic car- 
bon stars (Jura 1986), the range of gas mass-loss rates covered by 
the grid is 3.0 x 10" 10 M yr 1 to 4.3 x 10~ 5 M yr" 1 . We plot the 
GRAMS dust mass-loss rates against their luminosities in Fig. [2] 
The luminosities and mass-loss rate s determined for some LMC 
carbon stars that were m odeled by Ivan Loon et alJ d 19991) and 
IGroenewegen et all (120091) are also shown for comparison. The 
GRAMS grid has good coverage of the range of luminosities and 
mass-loss rates calculated by these studies, except at very high 
luminosities and mass-loss rates - in particu lar, there are three 
stars modeled by IGroenewegen et al.l J2009) with higher mass- 
loss rates. We would like to point out that this is no longer an 
issue if we adopt the same amorphous carbon dust opacities as 
in the Groenewegen et al. study (see Sect. l5.1.2l for details). The 
current optical depth limit of Tu.3 =4 was chosen based on the 
reddest sources in the SAGE study - with the current set of mod- 
els, we are able to produce redder mid-IR colors than those ob- 
served in the SAGE sample (see Sect. 14. 2\ . The choice was also 
based on results from our empirical study (Paper I). As mod- 
els with higher optical depth are computationally intensive, we 
decided to add optically thicker models in the future if required. 

4.2. Color-magnitude and color-color diagrams 

The figures in this section compare the GRAMS synthetic pho- 
tometry with SAGE O-rich, C-rich and extreme AGB candi- 
dates. These plots demonstrate that the GRAMS grid is able to 
reproduce the range of observed colors for carbon star candi- 
dates. We also include sources from the SAGE-Spec study that 
display m olecular features and/or dust signatures typical of car- 
bon stars dKemper et alJ2 010: Wo ods et alJ201 ll) . Only the sub- 
set of SAGE-Spec sources that had a full 5-37 urn spe ctrum are 
shown here. We refer the reader to iBlum etalJ (120061) for a de- 
tailed description of the various stellar populations observed in 
these diagrams and to Paper IV for discussion on the locations 
of the oxygen-rich AGB candidates and models. 

We would like to emphasize that comparisons on CMDs and 
CCDs are only shown in order to demonstrate the coverage of 
the observed colors and magnitudes, and we are careful not to 
over-interpret the agreement. Some of our models may be unre- 
alistic, yet others may correspond to unphysical parameter val- 
ues. While we have made an effort to avoid unrealistic models by 
filtering out those with very high dust temperatures, our aim was 
to span a large parameter space and good agreement with data 
does not necessarily validate a model in the grid. Similarly, good 
agreement with data on a color-color diagram does not necessar- 
ily imply that a model is a good fit to the data. Two models with 
significantly different values for the same parameter (e.g., lumi- 
nosity) but similar SED shapes might also end up close to each 
other on a color-color diagram, resulting in a large uncertainty 
for the parameter calculated for data based on just this diagram. 

4.2.1. K s vs. J-K s CMD 

This CMD is used to distinguish between stars with oxygen- 
rich and carbon-rich dust chemistries (see Paper I and references 
therein) because of a clear division of the stellar distribution into 
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Fig. 2. Top: The distribution of luminosities of LMC C-rich (solid) and extreme (dotted) AGB candidates from SAGE. The classical 
AGB luminosity limit is also shown (dashed line). Bottom: The range of luminosities and dust mass-loss rates covered by the 
GRAMS carbon star models. As discussed in Sect. 13.31 only models with T[ n < 1800 K ar e shown. Gray symbol s show the 
lumi nosities and mass-loss ra tes of LMC carbon stars as calculated from detailed modeling by Ivan Loon et al.l (1 19991) (triangles) 
and lGroenewegen et al.l (12009) (diamonds). 



two branches - the bluer branch, consisting of the O-rich AGB 
candidates, extends to brighter K s magnitudes where the red su- 
pergiants are located. The red branch reaches a maximum bright- 
ness with increasing color, at which point the near-infrared flux 
becomes progressively attenuated by the increasing amount of 
circumstellar dust. This branch consists of the C-rich and ex- 
treme AGB candidates. The carbon stars along the boundary be- 
tween the two branches {i.e., the bluest carbon stars) contain lit- 
tle or no dust around them; t he GRAMS models for low optical 
depths (as well as the bare Arin ger et al.l (|2009) photospheres) 
should lie along this boundary. Fig. [3] shows the K s vs. J-K s 
CMD for SAGE AGB candidates, with the model grid points 
superimposed. As expected, most of the A09 photospheres lie 
along the O/C boundary on this diagram. The NIR fluxes of AGB 
st ars are significant l y affe cted by pulsation; for instance, the data 
of IWhitelock et al] (120031) for selected LMC carbon stars sug- 
gests a median peak-to-peak variability of ~ 1 .0 mag in the J and 
K s bands. The agreement between the data and models is very 
good, considering that the A09 models do not include the effects 
of dynamical processes. A small number of A09 models possess 
much bluer colors and run across the O-rich AGB population. 
These are the warmer (T^ > 3500 K) photospheres which, as 
mentioned in Sect. l2.1.Tl are probably more representative of the 
post-AGB phase. Overall, the GRAMS grid provides excellent 
coverage of the J-K s colors observed for C-rich candidates and 



the extreme AGB candidates with detections in the near-infrared 
as well as carbon stars identified in the SAGE-Spec program. 
Photometry for OGLE LMC LPV 28579, which was studied in 
Paper III, is also shown. A series of vertical bands are visible in 
the models; these bands correspond to models with similar op- 
tical depths. For instance, the band with th j = 0.3 is seen just 
blueward of the photometry for OGLE LMC LPV 28579. 

lAringer et alj ((2009) noted that the H-K s colors predicted by 
their models were systematically bluer than observations. They 
explained the disagreement at warmer temperatures as probably 
arising from the scaling of their C2 opacity. We observe a sim- 
ilar disagreement between our models and SAGE data in the H 
band. For this reason, we do not discuss the 2MASS color-color 
diagram or CMDs involving the H band. 

4.2.2. [3.6] vs. J-[3.6] CMD 

This CMD, shown in Fig. @J is used in iBlum et aT] (120061) and 
Paper I to select extreme AGB stars. It is somewhat similar to 
the K s vs. J-K s CMD. The 3.6 fim magnitude is less affected by 
dust reddening and is therefore a better luminosity proxy for the 
extreme AGB candidates than the K s magnitude. The GRAMS 
grid shows good coverage of the entire range of J-[3.6] color, 
reproducing the colors of almost all the extreme AGB stars with 
J-band detections. Once again, the models form vertical bands 
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Fig. 3. The K s vs. J-K s CMD shows the model grid (light blue points) overlaid onto the SAGE AGB candidates from Paper I (gray: 
O-rich, orange: C-rich and pink: extreme). The yellow diamonds are SAGE sou rces identified as carb on-rich from their S AGE-Spec 
spectra. The large red triangle is OGLE LMC LPV 28579. Also shown are the Aring eret al.l ([2009) photospheres (green squares) 
used to generate the grid. The "spike" running across the O-rich AGB candidates consists of models with 7^ warmer than 3500 K. 
The cloud of models with Tn.3 = 0.3 is seen at J-K s ~ 3 mag. 




Fig. 4. The [3.6] vs. J-[3.6] CMD, same symbols as in Fig. [3] The cloud of models with Tn 3 = 0.3 is seen at J-[3.6] ~ 5 mag. 



of increasing optical depth. OGLE LMC LPV 28579 is seen 
to lie on the m.3 = 0.3 group of models. However, there are 
many C-rich candidates with J— [3.6] colors up to ~ 0.3 mag 
bluer than the model photospheres. This discrepancy in color 
can arise from misclassification of O-rich AGB stars as carbon- 
rich (resulting from our somewhat artificial definition of the O/C 
boundary based on near-IR colors) as well as the photometric un- 



certainty associated with the SAGE 3.6 fim photometry for our 
list of AGB candidates, which can be as high as 0.2-0.3 ma^3 



Also see the SAGE Data Products Description file available at 



http : //irsa . ipac . caltech . edu/data/SPITZER/SAGE/doc/SAGEDataPr 
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Fig. 5. The [24] vs. [8.0]-[24] CMD, same symbols as in Fig. [3] 



4.2.3. [24] vs. [8.0]-[24] CMD 

Fig. shows the [24] vs. [8.0]-[24] CMD. In Paper I, we used 
this diagram to illustrate the "bright" and "faint" O-rich AGB 
populations observed in SAGE data. The carbon stars are su- 
perimposed over most of the bright O-rich population, while 
the faint population is clearly seen at [8.0]-[24] <: 1 mag and 
[24] ^ 8 mag. The GRAMS models cover the entire range of 
MIPS 24 fim magnitudes observed for LMC carbon-stars. We 
note that there are several SAGE C-rich candidates as well as 
some SAGE-Spec carbon stars with [8.0]-[24] colors up to ~0.5 
mag bluer than the bluest GRAMS models (in fact, bluer than 
the A09 photospheres). While this discrepancy might potentially 
hint at limitations of the underlying stellar models used to gen- 
erate the grid, it is not inconsistent with the spread expected due 
to photometric uncertainties and variability. Using the uncertain- 
ties from Paper I for our C-rich and extreme AGB candidates, 
we estimate that the 3cr uncertainty in the [8.0]-[24] color can 
be up to about 0.55 mag (For more details, we refer the reader 
to the SAGE Data Products Description file). We calculate the 
median change in the color between th e two epochs of SAGE 
observations for the variable stars from IViih etal.l(l2009i) to be 
around 0. 1 mag. The observed discrepancy could thus be a com- 
bined effect of these uncertainties. 

There are also about 40 C-rich and extreme AGB candi- 
dates with [8]-[24] colors redder than the reddest models avail- 
able (~3.3 mag). We note here that at redder colors the ex- 
treme AGB candidate sample may be contaminated by young 
stellar objects (see discussion in Paper I). However, three of 
the very red sources mentioned above - IRAS 04535-6616, 
2MASS J0503 1662-6549450 and OGLE LMC LPV 16169 (not 
seen in Fig. [5] due to its extremely red color of 6.02 mag) - have 
been identified as carbon stars from their SAG E-Spec spectra. 
The fir st two are known C-rich AGB candidates: lvan Loon et alJ 
(119971) discovered the near-IR counterpart of IRAS 04535- 
6616 and classified it as a carbon star based on its colors, 
while Kontiza s~et alJ d200lh identified the carbonaceous nature 



of 2MASS J0503 1662-6549450 based on the detection of Swan 
C 2 bands in its optical spectrum. OGLE LMC LPV 16169 is 
seen projected against NGC 1835, a low-metallicity ([Fe/H] = 
-1.8) globular cluster with a considerable number of RR Lyrae 
variables (e.g., Walker 1993). As an extremely dust-enshrouded 
star, its presence in the metal-poor, low-mass population of 
NGC 1835 makes OGLE LMC LPV 16169 a very interesting 
object worth studying in detail. To be able to model such ex- 
tremely reddened sources, we will include higher optical depth 
models in future versions of the grid, and also consider more 
dust species. 

4.2.4. K s -[3.6] vs. J-K s CCD 

Fig. [6] shows the K s -[3.6] vs J-K s two-color diagram, which 
compares the emission from the warmest regions of the dust 
shell to that from the stellar photosphere. We show the GRAMS 
carbon-star grid superimposed on the SAGE AGB candidates 
from Paper I and the SAGE-Spec carbon stars on this diagram. 
The GRAMS models separate into groups of equal optical depth 
(in Fig. OGLE LMC LPV 28579 is seen to fall close to the 
group of models with optical depth thj = 0.3). This is more ob- 
vious at redder colors because of the coarse sampling at higher 
optical depths in the current version of the model grid. The J- 
K s color shows an interesting trend with effective temperature 
of the central star: for each set of models with the same opti- 
cal depth and inner shell radius, decreasing T e $ first causes a 
reddening in J-K s , due to the increasing continuum flux in the 
K s band relative to the J band. However, a lower value also 
means increased strength in the ~2.3 yum CO absorption feature, 
which eventually counteracts the increase in continuum K s flux 
and shifts the coolest models towards bluer J-K s . This trend with 
r e ff is depicted in Fig.|6]for the GRAMS models with logg = 0, 
M - 2 M , C/O = 2, T11.3 = 0.6 and R in = 7 R stm . At low optical 
depths (r < 0.1), the K s -[3.6] color is affected by the K s band 
CO absorption feature as well as strong HCN + C2H2 absorption 
near 3 jjm. With decreasing T^, the K s -[3.6] color first moves to 
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Fig. 6. The K s -[3.6] vs. J-K s CCD, same symbols as in Fig. [3] The "arc" slightly blueward of OGLE LMC LPV 28579 consists of 
the Tn.3 = 0.3 models. The variation of the J-K s color with effective temperature of the central star is depicted by the solid black 
curve for logg = 0, M = 2 M , C/O = 2, Tn 3 = 0.6 and R[ n = 7 R stal (T e ff decreases in the direction of increasing K s -[3.6] color). 




Fig. 7. The [3.6]-[4.5] vs. [5.8]-[8.0] CCD, same symbols as in Fig.H 



bluer values, reaches a minimum then reverses this trend. This 
turnover trend is not seen in either color for GRAMS models 
with the highest optical depths (r > 1) due to the large amounts 
of dust which overwhelms emission from the central star. For 
these models, a decrease in effective temperature is accompa- 
nied by a reddening in both colors. 

The GRAMS models are able to reproduce the range of ob- 
served colors in this diagram. The colors of the bare photo- 
spheres are similar to the bluest carbon-star candidates on the 



boundary with the O-rich candidates, as would be expected if 
these represented carbon stars with little or no dust. However, 
at redder J-K s colors, a significant fraction of C-rich and ex- 
treme AGB candidates are not covered by the models; in fact 
these sources are up to 0.8 mag bluer in K s -[3.6] than the mod- 
els. We have already mentioned that the photometric uncertainty 
in the 3.6 /mi SAGE data can be up to a few tenths of a mag- 
nitude. This uncertainty alone is not sufficient to explain the 
observed discrepancy. LMC carbon stars are known to exhibit 
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strong C2H2 absorption in the 3 pm region; as discussed previ- 
ously, this can depress the flux in the 3.6 pm band. At higher dust 
optical depths, emission from the circumstellar du st may some- 
what reduce this effect (see, e. e. . Ivan Loon et a which 
would explain why our high T\\ .3 models reproduce the colors 
of the reddest extreme AGB candidates on this plot. 

4.2.5. [5.8]-[8.0] vs. [3.6]-[4.5] CCD 

Fig. Q shows the IRAC-only two-color diagram. An increase in 
[3.6]-[4.5] color is accompanied by an increase in the [5.8]- 
[8.0] color for the extreme AGB candidates. This is simply a re- 
sult of increasing dust emission in the mid-IR due to increasing 
optical depth. A significant fraction of the C-rich AGB candi- 
dates, however, show a decreasing [5.8]-[8.0] color with redder 
[3.6]-[4.5] colors. This trend is also seen to a smaller extent in 
the distribution of the bluest ([3.6]-[4.5] < 0) model grid points. 
There is a strong CO absorption feature in the 5.8 pm band (see, 
e.g., Fig. 2 in Paper III) which could be filled in by emission 
from dust with increasing optical depth. 

The model colors agree well with those observed for the data, 
especially for the extreme AGB candidates. The photometric un- 
certainties as well as pulsation amplitudes in the IRAC bands are 
lower than in the near-IR, which may explain the smaller spread 
in the data around the models. At the bluest colors, there are 
many C-rich AGB candidates that are bluer than the model pho- 
tospheres. This could be explained, once again, on the basis of 
the absorption features in the IRAC 3.6 (HCN + C 2 H 2 ) and 5.8 
pm (CO) bands. 

We also note that the carbon-star models have continuous 
coverage of the extreme AGB candidates over the entire range of 
colors. Compare this to the oxygen-rich models from Paper IV, 
which do not reproduce the colors of the stars in the range 0.2 < 
[3.6]-[4.5] < 1 mag (see Fig. 5 and the discussion in Sect. 4.1.2 
in Paper IV). This difference in relative coverage can be used 
to identify the chemical types of the extreme AGB candidates 
through SED fitting. 

5. Discussion 

5.1. Preliminary SED fits 

5.1 .1 . Fits to the Ivan Loon et al.l ||1999|) sources 

We are ultimately interested in fitting the SEDs of the entire 
SAGE AGB star candidate list. In this section, we evaluate sim- 
ple chi-squared fits to the SEDs of spectroscopically confirmed 
carbon stars in the LMC for which detailed models already exist 
in the literature . Speci fical ly, we consider the sources studied by 
Ivan Loon et al.l (Il999l) and lGroenewegen et alJ J2009) (hereafter, 
vL99 and G09 respectively). We also compare the results of such 
SED fitting of the carbon star OGLE LMC LPV 28579 with 
the parameters derived from detailed modeling of this source in 
Paper III. The chi-squared fitting considered in this section al- 
lows a simple consistency check of the GRAMS grid. This tech- 
nique is able to predict accurate luminosities and dust mass-loss 
rates (see Paper III), but it may not be able to produce strong 
constraints on other parameters. For instance, we expect that a 
set of reasonable fits to a reddened source (e.g., a carbon star 
from the vL99 or G09 sample) will exhibit considerable degen- 
eracy in predicted stellar parameters. In what follows, therefore, 
we only compare the luminosities and dust mass-loss rates for 
the vL99 and G09 sources with the corresponding GRAMS best- 
fit values. We use the gas:dust ratio provided in vL99 and G09 



(500 and 200 respectively) to convert their total mass-loss rates 
to dust rates. Our u exp value is kept fixed at the same value as 
that of G09. However, vL99 use a luminosity-dependent expan- 
sion velocity given by 



lOkms- 1 \3xlO 4 L / 

For each vL99 source, we compare both the "scaled" (y exp as in 
Eq. [TJ and "unsealed" (i/ exp = 10 km s ) versions of the dust 
mass-loss rate with our best-fit value in Sect. 15.1.11 

Evolved AGB stars have significant variability amplitudes at 
visible wavelengths; this variation must be taken into account in 
order to reconstruct the SED from multi-wavelength photometry 
obtained at different pulsation phases. We will consider these de- 
tails in future papers. Presently we fit only the JHK S and Spitzer 
bands for these sources. As discussed in Sect. 14.11 our models 
provide substantial coverage of the range of observed luminosi- 
ties. Moreover, at the moment we would like to demonstrate the 
general agreement between fits from our model grid and those 
from detailed studies. For these reasons, we do not vary the 
model luminosities in order to find the best fit. We will treat the 
luminosity scale as an extra free parameter when we consider 
SED fitting in detail. 

5.1 .2. Fits to the Groenewegen et al. (2009) sources 

In their paper, van Loon et al. modeled 31 carbon stars (57 
LMC evolved stars total) with ISO photometry supplemented by 
ground-based near-IR photometry. ISO spectroscopy was also 
available for 10 of these stars. In Paper I, we used the vL99 mass- 
loss rates to estimate the total dust injection rate into the LMC 
from the entire AGB population. We found nearest-neighbor 
matches in our AGB candidate list (see Paper I) for all of the 
van Loon et al. carbon star£3- We fit GRAMS models to each of 
these stars, fitting only the near- and mid-IR bands as previously 
mentioned. Fig.[8]shows best-fit GRAMS spectra for four of the 
vL99 stars. Overall, our best-fit spectra reproduce the observed 
photometry well. The best-fit luminosities and dust mass-loss 
rates also show a good agreement with the values calculated by 
van Loon et al., as shown in Figs. [10] and QT| These plots demon- 
strate the power of the simple model grid fitting technique to 
obtain reliable luminosities and mass-loss rates. The GRAMS 
luminosities are within a factor of ~ 1.75 of the vL99 values. 
The median discrepancy is -13%, with a standard deviation of 
29% around this value. With the exception of TRM 45, for which 
we obtain a poor fit, our mass-loss rates agree with those of vL99 
(unsealed) to within a factor of 4.5. Our rates are discrepant from 
the unsealed vL99 rates by a median value of 8% with a standard 
deviation of 83% and from the scaled versions by -13% with a 
standard deviation of 125%. 

Groenewegen et al. used Spitzer IRS spectra to model 101 
Magellanic Cloud carbon stars, using multi-epoch photometry 
to constrain the optical variability. Their dataset included 68 
LMC carbon stars. We found counterparts in our list of AGB 
candidates for all but three of the G09 stars - MSX LMC 95, 
MSX LMC 1384 and NGC 1978 IR1. However, on searching the 
SAGE Archive point source lists, we were able to find a match 
for MSX LMC 95 in the IRAC Epoch 1 Archive. We also found 
detections for MSX LMC 1384 and NGC 1978 IR1 in the IRAC 



14 We excluded IRAS 05289-6617 from further analysis because it 
lacked detections in the optical an d near-IR bands as w ell as most of 
the Spitzer bands (see discussion in lvan Loon etal.lll999h . 
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Fig. 8. GRAMS fits to the SEDs of four carbon stars modeled bv lvan Loon et all {l999). 



Mosaic Photometry catalog. Fig. [9] shows t he best-fit GRAMS 
spectr a for four carbon stars studied by iGroenewegen et ail 
(2009). The best-fit luminosities and dust mass-loss rates for the 
entire set are compared with the G09 estimates in Figs.lTOlandlTTl 
respectively. Once again, the luminosities predicted by GRAMS 
SED fitting are within a factor of ~ 1 .75 of the values determined 
from detailed modeling. The median discrepancy is very low at 
-1%, with a standard deviation of 31%. However, the behavior 
of the mass-loss rates is different from that observed with the 
vL99 sample. 

While the luminosities are more or less evenly distributed 
about the 1:1 line for GRAMS fits to both the vL99 and G09 
sources, our dust mass-loss rate estimates for most of the G09 
sources are systematically lower than theirs by a factor of about 
2—4. The median discrepancy in this case is -168%, with a stan- 
dard deviation of almost 200% around this value. This discrep- 
ancy is due to our different choice of am orphous carbon optical 
constants. The lRouleau & Martini d 19911) optical constants used 
in G09 give rise to a consistently lower absor ption effic iency 
over the ~0. 1-100 fim range (see, e.g., Fig. 1 of Suh 2000), thus 
requiring a higher dust shell mass in order to fit a given mid- 
IR flux. This in turn increases the calculated mass-loss rate. In 
order to verify the effect of this change, we performed detailed 
radiative transfer calculations on ten of the G09 sources using 
the GRAMS best-fit paramete rs from SED fitting as input, with 
two modifications: we used the Ro uleau & Martini (I 1 99 ll) optical 
constants for amorphous carbon and set the SiC mass fraction to 
the value determined in G09 for each star. We found that in gen- 



eral the iRouleau & Martini (1199 ll) optical constants gave mass- 
loss rates that were 2-4 times higher than those predicted by the 
GRAMS best fits, explainin g the systematic offset observed in 
Fig-El For comparison, the Preibisch et alj dl99 3) amorphous 
carbon optical constants used bv lvan Loon et al.l ( 1999) ha ve ab- 
sorption efficiencies similar to that of the Zubko et al. (1996) 
dust model. We computed 2Dust mo dels for five of the vL99 
sources with the Preibisc h et al.l (1 19931) dust constants and found 
that the resulting dust mass-loss rates were at most 12% differ- 
ent from those predicted by the corresponding GRAMS best-fit 
models. In developing this first version of the model grid, we 
used a single set of dust properties. We will investigate the effect 
of varying the dust properties in detail in future versions of the 
grid. 

TRM 88 was modeled by both vL99 and G09. Fig.fl2lshows 
the GRAMS best-fit spectrum to the observed photometry of 
TRM 88. The GRAMS best-fit luminosity ( 1 1 700 L Q ) and mass- 
loss rate (3.4 x 10~ 9 M yr 1 ) are closer to the values estimated 
by van Loon et al. In this case, the GRAMS fit is somewhat poor 
as it under-predicts the mid-IR emission resulting in a low mass- 
loss rate. The marginal quality of the fit is partly due to the fact 
that we have not taken the near-IR variability into account; the 
fluxes are currently weighted by their photometric errors which 
are quite small. The best-fit model has an optical depth T11.3 = 
0. 1 and Rja = 12 R stm . We do not have enough resolution in opti- 
cal depth (the closest higher optical depth is 0.2) and inner radius 
(we only explore four values) to provide a better fit to the data. 



13 



S. Srinivasan et al.: The GRAMS Carbon-Star Model Grid 





Fig. 9. GRAMS fits to the SEDs of four carbon stars modeled bv lGroenewegen et all (|2009). 




Fig. 10. Luminosities calculated from GRAMS SED fitting plot- 
ted against the l umin osities obtained from detailed mode l ing by 
Ivan Loon et al.l d 19991) (squares) and lGroenewegen et al.1 12009) 
(circles). The solid line represents a 1:1 agreement. 



Fig. 11. Same as Fig. \10\ but for dust mass-loss rates. Also 
shown are lines along which the GRAMS mass-loss rate is lower 
than the rates from detailed modeling (vL99 or G09) by a factor 
of 2 and 4 (dotted and dashed line respectively). 
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Fig. 12. GRAMS best-fit for TRM 88. The best-fit luminosity 
and mass-loss rates are more consistent with the vL99 estimates. 



Despite these details, we find good agreement with the vL99 re- 
sults. 

5.1 .3. Fit to OGLE LMC LPV 28579 
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Fig. 13. Best fit SEDs to the SAGE photometry (dots) and 
SAGE-Spec IRS spectrum for OGLE LMC LPV 28579. The 
solid curve is a result of the detailed model from Paper III, while 
the dashed curve is the GRAMS best-fit. A comparison of the fit 
parameters is also shown. 

In Paper III we presented a 2Dust model best-fit model for 
OGLE LMC LPV 28579. The dust properties determined for this 
best-fit model were subsequently used for the entire carbon-star 
grid. Our detailed model suggested a luminosity of 6580 L for 
the SAGE Epoch 1 observations, an optical depth of T11.3 = 0.27, 



and a dust mass-loss rate of 2.5 x 10 9 M yr The mass-loss 
rate derived for OGLE LMC LPV 28579 was found to be con- 
sistent with the rates derived from various empirical relation- 
ships (see discussion in Paper III for details). In this paper we 
fit the observed SED for OGLE LMC LPV 28579 using models 
from the GRAMS grid. This serves as a consistency check for 
the model grid. The best-fit GRAMS SED gives a luminosity of 
6170 L , T11 3 = 0.3 and a dust mass-loss rate of 2.4 x 10~ 9 M 
yr -1 . These numbers are within a few percent of the values cal- 
culated in Paper III. Fig. [13] compares the best-fit models from 
both papers with the SAGE Epoch 1 photometry and SAGE- 
Spec spectrum for OGLE LMC LPV 28579. We also obtain good 
agreement with the values obtained for the optical depth, dust 
shell inner radius as well as temperature at the inner radius. 

6. Summary 

We constructed a grid of carbon star models using the radiative 
transfer code 2Dust and for a range of various stellar and dust 
shell parameters. We intend to use these models in conjunction 
with our models for O-rich AGB and red supergiant stars de- 
scribed in Sar gent et al.l d201 ll) to investigate the mass-loss re- 
turn from evolved stars to the LMC. The models can also be 
used for similar estimates from large photometric samples of C- 
rich AGB stars. The grid covers luminosities from 2000 L to 
26 000 L by using model photospheres spanning temperatures 
in the range 2600^1-000 K. Assuming spherically symmetric dust 
shells, and constant expansion velocity and mass-loss rate, we 
perform radiative transfer using the 2Dust code for 1 1.3 //m op- 
tical depths ranging from 10~ 3 to 4 and values for the inner radii 
of 1.5, 3, 4.5, 7 and 12 times the stellar radius. This results in 
over 12 000 models with dust temperatures under 1800 K, with 
dust mass-loss rates in the range 10~ 12 - 10~ 7 M yr -1 . 

We synthesize photometry for these models in optical as well 
as near- and mid-infrared bands. The entire set of models, in- 
cluding spectra and synthetic photometry, will soon be available 
at the 2Dust website at STScI. We compare the resulting colors 
and magnitudes with those observed for AGB candidates in the 
SAGE survey and confirmed AGB stars from the SAGE-Spec 
survey, finding good overall agreement with these data. Using 
a chi-squared fitting routine, we obtain best-fit spectra, lumi- 
nosities and ma ss-loss rates for s pectro scop ically identified car- 
bon st ars in the Ivan Loon et alJ (Q999) and iGroenewegen et ail 
(l2009h samples. The luminosities predicted from simple fitting 
to the photometry are in good agreement with those determined 
by these detailed models. Our mass-loss rates for the van Loon 
et al. sample agree well with their values. However, our rates 
for the Groenewegen et al. sample are lower by a factor of 2-4, 
most likely due to a different choice of dust properties. We find 
excellent agreement between the results of detailed modeling of 
OGLE LMC LPV 28579 from Paper III and the fitting employed 
in this work. 

One of the aims of our study is to provide a general-use fitter 
that can be applied to large sets of photometric data. This will be 
the focus of one of our future papers. The GRAMS grid will en- 
able the assessment of mass-loss return from galaxy-wide point 
source catalogs from projects such as SAGE. While most ex- 
treme AGB stars are probably carbon-rich, some are very bright 
OH/IR stars. Currently, there is no way to clearly distinguish 
between these sources with photometry alone. We hope that a 
comparison of our results with ongoing studies such as AKARI 
and WISE that include filters sensitive to silicate and SiC fea- 
tures will enable us to define the separation between C-rich and 
O-rich sources in this extreme regime. 
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